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Abstract 

Three series of tensile relaxation tests are performed on isotactic polypropylene 
at room temperature in the vicinity of the yield point. In the first series of exper- 
iments, injection-molded samples are used without thermal pre-treatment. In the 
second and third series, prior to testing the specimens are annealed at 130 °C for 4 
and 24 hours, respectively. 

Constitutive equations are derived for the time-dependent response of semicrys- 
talline polymers at isothermal loading with small strains. A polymer is treated as 
an equivalent temporary network of macromolecules bridged by junctions (physical 
cross-links, entanglements and crystalline lamellae). Under loading, junctions slip 
with respect to their positions in the bulk material (which reflects the viscoplas- 
tic behavior), whereas chains separate from their junctions and merge with new 
ones at random times (which reflects the viscoelastic response). The network is 
thought of as an ensemble of meso-regions (MR) with various activation energies 
for detachment of chains from temporary nodes. 

Adjustable parameters in the stress-strain relations are found by fitting obser- 
vations. Experimental data demonstrate that the shape of the relaxation spectrum 
(characterized by the distribution of MRs with various potential energies) is indepen- 
dent of mechanical factors, but is altered at annealing. For specimens not subjected 
to thermal treatment, the growth of longitudinal strain does not affect the volume 
fraction of active MRs and the attempt rate for detachment of chains from their 
junctons. For annealed samples, the concentration of active MRs increases and the 
attempt rate decreases with strain. These changes in the time-dependent response 
are attributed to broadening of the distribution of strengths of crystalline lamellae 
at annealing. 

1 Introduction 

This paper is concerned with the influence of annealing at an elevated temperate on the 
nonlinear viscoelastic response of isotactic polypropylene (iPP) at room temperature. The 
objective of this study is three-fold: 
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1. to report experimental data in tensile relaxation tests on specimens annealed for 
various amounts of time at strains in the vicinity of the yield point, 



2. to derive stress-strain relations for the nonlinear viscoelastic behavior of a semicrys- 
talline polymer at isothermal uniaxial deformation, 

3. to assess the effect of annealing on the time-dependent response of iPP in terms of 
the constitutive model. 

Isotactic polypropylene is chosen for the analysis because of numerous applications of 
this polymer in industry (oriented films for packaging, reinforcing fibres, nonwoven fab- 
rics, blends with thermoplastic elastomers, etc.). The goal of this study is to establish 
some correlations between mechanical properties, morphology and processing conditions 
(annealing at an elevated temperature) for injection- molded specimens. For a review of 
previous works on this subject, see JH and the bibliography therein. 

The nonlinear viscoelastic response of polypropylene was analyzed by Ward and Wolfe 
0], see also 0, and Smart and Williams |Q three decades ago, and, more recently, by 
Ariyama ^, |^, |^, Wortmann and Schulz |TU|, Ibhadon |TI|, Tomlins |]T^, Read and 



Tomlins 
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p]5|] , Dutta and Edward [|I^ and Tomlins and Read 
The effect of physical aging (anneahng at an elevated temperature followed by quench 
to ambient temperature) on the time-dependent behavior of PP was studied by Struik 
17| , |T^ , Chai and McCrum |jl9|, Ibhadon [|ll|], Tomlins [jl2|. Read and Tomlins 0, |15 
and Tomlins and Read . 



Dynamic mechanical analysis shows that the loss tangent of iPP demonstrates two 
pronounced maxima being plotted versus temperature The first maximum (/3- 

transition in the interval between T = —20 and T = 10 °C) is associated with the glass 
transition in the most mobile part of the amorphous phase, whereas the other maximum 
(a-transition in the interval between T = 70 and T = 110 °C) is attributed to the glass 
transition in the remaining part of the amorphous phase. This conclusion is confirmed by 
DSC (differential scanning calorimetry) traces for quenched PP that reveal an endoterm 
at T = 70 °C which can be ascribed to thermal activation of amorphous regions with 



restricted mobility under heating |21 



Isotactic polypropylene exhibits three different crystallographic forms: monoclinic a 
crystallites, (pseudo) hexagonal (3 structures, orthorhombic 7 polymorphs, and "smectic" 
mesophase (arrays of chains with a better order in the longitudinal than in transverse 
chain direction). For a detailed review of iPP polymorphs, the reader is referred to 
the survey [^]. At rapid cooling of the melt (which is typical of injection molding), a 
crystallites and smectic mesophase are mainly developed, whereas metastable jS and 7 
structures arise as minority components. Crystallization of (3 forms occurs either under 
stresses or with the help of /3-nucleating agents added to the melt 



I, P4|. Formation of 7 

polymorph requires high pressure for commercial grades of iPP, while it can observed at 
atmospheric pressure in isotactic polypropylene with low molecular weight pSf . A unique 
feature of a structures in iPP is the lamellar crosshatching: development of transverse 
lamellae in spherulites that are oriented in the direction perpendicular to the direction of 
radial lamellae . 
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Scanning electron microscopy |2^, atomic force microscopy |2S[ and X-ray diffrac- 
tion ||21], ^|2^ reveal that in injection-molded specimens a spherulites have the character- 
istic size of the order of 100 /xm and they contain crystalline lamellae with thickness of 10 
to 20 nm. The amorphous phase is located between spheruhtes and inside the spherulites 
between lamellae. It consists of (i) relatively mobile chains between spherulites and be- 
tween radial lamellae inside spherulites, and (ii) severely restricted chains in the regions 
bounded by radial and tangential lamellae in a spherulites. 

Annealing of injection-molded iPP at an elevated temperature results in (i) secondary 
crystallization of a part of the amorphous phase, (ii) thickening of radial lamellae, (iii) 
development of subsidiary lamellae, (iv) formation of lamellar superstructure, and (v) 
growth of the crystal perfection ||30|. Other changes in the crystalline morphology of 
iPP driven by thermal treatment are the subject of debate. Some researchers [|^, 



28| , [29|] conclude that the fraction of (5 spherulites increases at annealing in the interval 
of temperatures between 110 and 140 °C, which enhances ductility of iPP and improves 
its impact properties. According to other authors ETI, EB, pT|, annealing of iPP induces 



transformation of the smectic phase into monoclinic ol spherulites without noticeable 
development of /5 polymorph. 

Mechanical loading results in inter-lamellar separation, rotation and twist of lamellae, 
fine and coarse slip of lamellar blocks and their fragmentation |2l[]. Straining of iPP 



specimens causes chain slip through the crystals, sliding and breakage of tie chains and 
activation of restricted amorphous regions driven by lamellar disintegration. In the post- 
yield region, these changes in the micro-structure imply cavitation, breakage of crystalls, 
and formation of fibrills p2 |. 



It is hard to believe that these morphological transformations in iPP can be adequately 
described by a constitutive model with a small number of adjustable parameters. To 
develop stress-strain relations, we apply a method of "homogenization of micro-structure" 
|33| . According to this approach, an equivalent phase is introduced whose deformation 



captures essential features of the response of a semicrystalline polymer with a complicated 
micro-structure. In this study, an amorphous phase is chosen as the equivalent phase 
because of the following reasons: 

1. The viscoelastic response of semicrystalline polymers is conventionally associated 
with rearrangement of chains in amorphous regions p9| . 



2. Sliding of tie chains along and their detachment from lamellae play the key role in 
the time-dependent response of semicrystalline polymers |24|, p5|] . 



The viscoplastic flow in semicrystalline polymers is assumed to be "initiated in the 
amorphous phase before transitioning into the crystalline phase" p6|. 



Conventional models for polyethylene |33|, polypropylene and poly(ethylene 



terephthalate) treat these polymers as networks of macromolecules. 

Above the glass transition temperature for the mobile amorphous phase, isotactic 
polypropylene is thought of as a network of chains bridged by junctions. Deformation of 
a specimen induces slip of junctions with respect to their positions in the bulk material. 
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Sliding of junctions reflects slippage of tie molecules along lamellae and fine slip of lamellar 
blocks which are associated with the viscoplastic behavior of a semicrystalline polymer. 



With reference to the concept of transient networks ^ ^ the viscoelas- 
tic response of iPP is modelled as separation of active chains from their junctions and 
attachment of dangling chains to temporary nodes. The network of macromolecules is 
assumed to be strongly inhomogeneous (this heterogeneity reflects the effect of spherulites 
on rearrangement of surrounding chains) , and it is treated as an ensemble of meso-regions 
(MR) with various potential energies for detachment of active strands. Two types of MRs 
are distinguished: (i) active domains where strands separate from junctions as they are 
thermally agitated (these MRs model a mobile part of the amorphous phase), and (ii) 
passive domains where detachment of chains from junctions is prevented. Passive MRs 
are associated with a part of the amorphous phase whose mobility is restricted by (i) 
radial and tangential lamellae and (ii) surrounding macromolecules (because of density 
fluctuations in the amorphous phase). 

Separation of active chains from temporary nodes is treated as a thermally-activated 
process whose rate obeys the Eyring equation with a strain-dependent attempt rate. 
An increase in the relaxation rate of amorphous polymers at straining is conventionally 
attributed to the mechanically-induced growth of "free volume" between macromolecules 
which, in turn, implies an increase in their mobility H, H, H. 

Stretching of a specimen results in (i) a mechanically-induced changes in the rate of 
detachmens of strands in active MRs and (ii) an increase in the concentration of active 
MRs. The latter is ascribed to (i) partial release of the amorphous phase in passive meso- 
domains driven by fragmentation of lamellae and (ii) breakage of van der Waals links 
between compactly packed chains in meso-domains with higher density. 

The exposition is organized as follows. Section 2 is concerned with the description of 
the experimental procedure. Kinetic equations for sliding of junctions and reformation of 
active strands are developed in Section 3. Constitutive equations for uniaxial deformation 
are derived in Section 4. In Section 5 these relations are applied to fit experimental data. 
A brief discussion of our findings is presented in Section 6. Some concluding remarks are 
formulated in Section 7. 



2 Experimental procedure 

Isotactic polypropylene (Novolen llOOL) was supplied by BASF (Targor). ASTM dumbell 
specimens were injection molded with length 14.8 mm, width 10 mm and height 3.8 mm. 
Three series of tests were performed. In the first series, the samples were used as received 
without thermal pre-treatment. In the second series, the specimens were annealed in an 
oven at the temperature 130 °C for 4 h and slowly cooled by air. In the third series of 
experiments, the specimens were annealed at the same temperature for 24 h and cooled 
by air. 

Differential scanning calorimetry measurements were carried out on STA 449/Netzsch 
apparatus at the heating rate 5 K/min. The specimens with weight of about 15 mg were 
tested in AI2O2 pans covered by lid. The thermal analyzer was calibrated with 7 references 
ranging from In to Ni. The specific enthalpy of melting, Aifm, equals 86.9, 98.1 and 101.7 
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J/g, for non- annealed specimens, specimens annealed for 4 h, and specimens annealed for 
24 h, respectively. With reference to we accept the value 209 J/g as the enthalpy of 
fusion for a fully crystalline polypropylene. The degree of crystallinity, Kc, is estimated 
as 41.6%, 46.8 % and 48.7 % for specimens not subjected to thermal treatment, samples 
annealed for 4 h, and specimens annealed for 24 h, respectively. 

Although the degree of crystallinity changes rather weakly (but consistently) with 
an increase in the annealing time, the shape of DSC curves is noticeably altered in the 
interval of temperatures between 120 and 160 °C. DSC traces depicted in Figure 1 are 
similar to those found by other authors, see, e.g., Labour et al. |23| attributed 



the growth of the low-temperature shoulder on the melting curve to the a — >■ /3 transition 



at annealing. According to lijima and Strobl |2g], the endothermal contribution on the 
low-temperature side of the DSC trace indicates the presence of a crystallites with varying 
stability. 

Uniaxial tensile relaxation tests were performed at room temperature on a testing 
machine Instron-5568 equipped with electro-mechanical sensors for the control of longi- 
tudinal strains in the active zone of samples (the distance between clips was about 50 
mm). The tensile force was measured by the standard loading cell. The engineering 
stress a was determined as the ratio of the axial force to the cross-sectional area of the 
specimens in the stress-free state. 

Any series of mechanical experiments included 9 relaxation tests at the longitudinal 
strains ei = 0.02, ea = 0.04, eg = 0.06, = 0.08, eg = 0.10, eg = 0.12, 67 = 0.14, eg = 0.16, 
69 = 0.18, which corresponded to the domain of nonlinear viscoelasticity, sub-yield and 
post-yield regions for isotactic polypropylene (the yield strain, ey, was estimated by the 
supplier as 0.13). Mechanical tests were carried out at least one day after annealing of 
specimens to avoid the influence of physical aging on the time-dependent response of iPP. 

Each relaxation test was performed on a new sample. No necking of specimens was 
observed in experiments (except for the test with eg = 0.18 on a specimen not subjected 
to thermal treatment, which was excluded from the consideration). In the kth relaxation 
test {k = 1, ... ,9), a specimen was loaded with the cross-head speed 5 mm/min (that 
roughly corresponded to the strain rate eo = 0.05 min~^) up to the longitudinal strain e^, 
which was preserved constant during the relaxation time = 20 min. 

The engineering stresses, a, at the beginning of the relaxation tests are plotted in 
Figure 2 together with the stress-strain curves for the specimens strained up to eg = 0.18. 
The figure demonstrates fair repeatability of experimental data. 

Figure 2 shows that annealing for 4 h results in a pronounced increase in stress com- 
pared to virgin specimens. The growth of the annealing time implies a decrease in stress 
in the sub-yield region. The discrepancy between the stress-strain curves for specimens 
annealed for 4 and 24 h practically disappears in the post-yield domain. 

Despite the coincidence of the stress-strain diagrams in the post-yield region for spec- 
imens annealed for 4 and 24 h, these samples demonstrate a noticeably different necking 
behavior. Necking of specimens not subjected to thermal pre-treatment occurs at the 
strain = 0.18, necking of specimens annealed for 4 h takes place at the strain = 0.25, 
whereas no necking is observed for specimens annealed for 24 h at stretching up to the 
strain e = 0.30. 

The longitudinal stress, a, is plotted versus the logarithm (log = log^g) of time t (the 
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initial instant t = corresponds to the beginning of the relaxation process) in Figures 3 
to 11. These figures demonstrate that the time of annealing strongly affects the shape of 
relaxation curves (especially, in the sub- yield region, see Figures 3 and 4). For any strain 
e > 0.02, stresses in the annealed specimens exceed those in the samples not subjected 
to thermal pre-treatment. In the sub-yield domain (e < 0.1) stresses in the specimens 
annealed for 4 h are higher than stresses in the samples annealed for 24 h. 

Our aim now is to develop constitutive equations for the time-dependent behavior of a 
semicrystalline polymer to be employed for the quantitative analysis of these observations. 



3 A micro- mechanical model 

A semicrystalline polymer is treated as a temporary network of chains bridged by junc- 
tions. The network is modelled as an ensemble of meso-regions with various strengths 
of interaction between macromolecules. Two types of meso-domains are distinguished: 
passive and active. In passive MRs, inter-chain interaction prevents detachment of chains 
from junctions, which implies that all nodes in these domains are permanent. In active 
MRs, active strands (whose ends are connected to contiguous junctions) separate from 
the temporary junctions at random times when they are thermally agitated. An active 
chain whose end detaches from a junction is transformed into a dangling chain. A dan- 
gling chain returns into the active state when its free end captures a nearby junction at 
a random instant. 

Denote by X the average number of active strands per unit mass of a polymer, by Xa 
the number of strands merged with the network in active MRs, and by Xp the number 
of strands connected to the network in passive MRs. Under stretching some crystalline 
lamellae (restricting mobility of chains in passive MRs) break, which results in a growth 
of the number of strands to be rearranged. As a consequence, the number of strands in 
active MRs increases and the number of strands in passive meso-domains decreases. This 
implies that the quantities Xa and Xp become functions of the current strain, e, that obey 
the conservation law 

Xa(e)+Xp(6)=X (1) 

Rearrangement of strands in active MRs is thought of as a thermally activated process. 
The rate of detachment of active strands from their junctions in a MR with potential 
energy u in the stress-free state of a specimen is given by the Eyring equation 



Facxp 



UJ 



where /cb is Boltzmann's constant, T is the absolute temperature, and the pre-factor Fa 
is independent of energy u and temperature T. Introducing the dimensionless activation 
energy u = oj/ik-oTo), where Tq is a reference temperature, and disregarding the effects 
of small increments of temperature, AT = T — Tq, on the rate of separation, F, we arrive 
at the formula 

F = Faexp(-a;). (2) 

We suppose that Eq. (2) remains valid for an arbitrary loading process, provided that 
the attempt rate. Fa, is a function of the current strain. Fa = Fa(e). 
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The distribution of active MRs with various potential energies is described by the 
probabihty density p{uj) that equals the ratio of the number, Na,{e,u), of active meso- 
domains with energy cu to the total number of active MRs, 

N,{e,uj) = X,{e)p{uj). (3) 

The distribution function for potential energies of active MRs, is assumed to be 

strain-independent. 

The ensemble of active mcso-domains is described by the function n,^{t^T,uo) that 
equals the number of active strands at time t (per unit mass) belonging to active MRs 
with potential energy uj that have last been rearranged before instant r e [0,t]. In 
particular, na(0, 0, uj) is the number (per unit mass) of active strands in active MRs with 
potential energy a; in a stress-free medium, 

na(0,0,cu) = Ara(0,a;), (4) 

and na(t, t, oo) is the number (per unit mass) of active strands in active MRs with potential 
energy uj in the deformed medium at time t (the initial time t = corresponds to the 
instant when external loads are applied to the polymer), 

n^{t,t,u)^N^{e{t),u). (5) 

The amount 



-(t,r,u;J 



dr 



Or 

equals the number (per unit mass) of dangling strands in active MRs with potential energy 
u! that merge with the network within the interval [r, r -|- dr] , and the quantity 



dr 



(t, r, uj)dT 



is the number of there strands that have not detached from temporary junctions during 
the interval [t, t]. The number (per unit mass) of strands in active MRs that separate 
(for the first time) from the network within the interval \t,t + dt] reads 

''-{t,0,u)dt, 



dt 

whereas the number (per unit mass) of strands in active MRs that merged with the 
network during the interval [r, r + dr] and, afterwards, separate from the network within 
the interval [t,t + dt] is given by 

-[t, T, U!)dtdT. 



dtdr 

The rate of detachment, F, equals the ratio of the number of active strands that separate 
from the network per unit time to the current number of active strands. Applying this 
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definition to active strands that merged with the network during the interval [r, r + dr] 
and separate from temporary junctions within the interval [t, t + dt], we find that 



QIQt^^^^^^^) ^ -r(e(t),a;) — (t,T,u;). 



(6) 



Changes in the function na{t,0,uj) are governed by two processes at the micro-level: 
(i) detachment of active strands from temporary nodes, and (ii) mechanically-induced 
activation of passive MRs. The kinetic equation for this function reads 



Oti d TV da 

-Q^{t,0,u;) = -r{e{t),u;)n,{t,0,u;) + -^{e{t),cu)-{t). 

The solution of Eq. (7) with initial condition (4) is given by 



na(t,0,a;) 



N^{0,uj) exp 
rt dN.. 



+ 



f 

Jo 5e ^ ^ ' ' dt 



T{e{s),u!)ds 



de 



r{e{s),uj)ds 



dr. 



It follows from Eq. (6) that 

dn 



dr 



{t, T, oj) = (p{t, oj) exp 



j T{e{s),uj)ds 



where 



dr 



{t,T,u) 



t=T 



To determine the function if{t,u!), we use the identity 

ft Qj^ 

n^{t,t,uj) = n^{t,0,uj) + -^{t,T,uj)dT. 
Equations (5) and (11) imply that 

ft Qrj^ 

na(i,0,a;) + -^{t,r,cu)dr ^ N^t),^)- 
Differentiating Eq. (12) with respect to time and using Eq. (10), we obtain 

^(t^^) + ^(^'0'^) + / ^(t,T,uj)dT = ^(e(t),^)^(t). 



dtdr ' ' de 
This equality together with Eqs. (6), (7) and (11) results in 

(p{t,u;) = r{e{t),u)n^{t,t,u;). 



(7) 



(8) 

(9) 
(10) 

(11) 
(12) 



(13) 



Substituting expression (13) into Eq. (9) and using Eq. (5), we arrive at the formula 



dr 



(i,T,a;) = r(e(i),a;)A^a(e(^),^)exp 



r(e(s), a;)ds 



(14) 



The kinetics of rearrangement of strands in active MRs is described by Eqs. (2), (3), (8) 
and (14). These relations are determined by (i) the distribution function p{iu) for active 
MRs with various potential energies co, (ii) the function ra(e) that characterizes the effect 
of strains on the attempt rate, and (iii) the function 

«a(6) = (15) 

that reflects mechanically-induced activation of passive MRs. 

Rearrangement of strands in active MRs reflects the viscoelastic response of a semicrys- 
talline polymer. The viscoplastic behavior is associated with the mechanically-induced 
slippage of junctions with respect to their positions in the bulk material. 

Denote by €a{t) the average strain induced by sliding of junctions between macro- 
molecules (the subscript index "u" means that e^{t) coincides with the residual strain 
in a specimen which is suddenly unloaded at instant t). The elastic strain (that reflects 
elongation of active strands in a network) is denoted by ee{t). The strains ee{t) and e^it) 
are connected with the macro-strain e{t) by the conventional formula 

eit)^e,it) + e^it). (16) 

We adopt the first order kinetics for slippage of junctions with respect to the bulk material, 
which implies that the increment of the viscoplastic strain, dsu, induced by the growth 
of the macro-strain, e, by an increment, de, is proportional to the absolute value of the 
stress (7, 

^ = 5Hsignf4), (17) 



de ^ dt 

where the pre-factor 5 is a non-negative function of stress, strain and the strain rate. 



de-. 



The last multiplier in Eq. (17) determines the direction of the viscoplastic flow of junc- 
tions. Equation (17) is presented in the form 

^{t) = B{a{t),e{t), ^)\a{t)\ sign[a{t)^(t)]^(t), 6,(0) = 0, (18) 



which implies that the rate of sliding vanishes when the macro-strain, e, remains constant. 



4 Constitutive equations 

An active strand is modelled as a linear elastic solid with the mechanical energy 

w{t) = ^i,e\t), 

where /i is the average rigidity per strand and e is the strain from the stress-free state to 
the deformed state. 
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For strands belonging to passive meso-domains, the strain e coincides with ee. Mul- 
tiplying the strain energy per strand by the number of strands in passive MRs, we find 
the mechanical energy of meso-domains where rearrangement of chains is prevented by 
surrounding lamellae, 

W,it) = i/xXp(6(t))62(t). (19) 



With reference to the conventional theory of temporary networks we assume that 
stresses in dangling strands totally relax before these strands merge with the network. 
This implies that the reference (stress-free) state of a strand that is attached to the 
network at time r coincides with the deformed state of the network at that instant. For 
active strands that have not been rearranged until time t, the strain e(t) coincides with 
ee(t), whereas for active strands that have last been merged with the network at time 
T E [0, t], the strain e{t, r) is given by 

e(t,r) = ee(t) - ee(r). 

Summing the mechanical energies of active strands belonging to active MRs with various 
potential energies, u, that were rearranged at various instants, r, we find the mechanical 
energy of active meso-domains. 



it,T,uj) ee(t) - ec(r) 



dr 



(20) 



The mechanical energy per unit mass of a polymer reads W(t) = Wg^^t) + Wp(t). Substi- 
tuting expressions (19) and (20) into this equality and using Eq. (16), we arrive at the 
formula 



W{t) 



(t,r,a;)((e(t)-e,(t) 



dr 



duj 



e r 



na(t,0,cu)(e(t) -eu(t) 



e,, r 



dr 



Differentiation of Eq. (21) with respect to time results in 

dW 



dt 



■(t)=/i 



(21) 



(22) 



where 



A{t) 



X,{e{t))[e{t) - e^{t) 
rt dn 



+ 



dr 



duj\n^(t,Q,uj) e{t) - e^{t) 

( f)r)^ 

duo 



(t,r,^)[(e(t)-e„(t))-(e(r)-e,(r)) 

dn 



dt 



(t,0,cu) t{t)-t^{t) 



+ 



dtdr 



(t,r,^)[(e(t)-e,(t)) -(e(r)-e,(r)) 



dr 



(23) 
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Bearing in mind Eqs. (1), (3), (5) and (11), we transform the first equality in Eq. (23) 
as follows: 



A{t) = X[e{t) - e^{t)\ - cLu -^{t,T,iu)[e{T) - e^ir) 



dr. 



(24) 



Substitution of expressions (1), (3), (6) and (7) into the second equality in Eq. (23) yields 

Jo 



+ 



V{e{t),u)duj{n^{t,Q,uj) e{t) - e^{t) 

* On 



dr 



'-{t,T,Lj)[[e{t)-e,{t)) - (e(T) -e,(T)) 



'dry 



(25) 



For uniaxial loading with small strains at the reference temperature Tq, the Clausius- 
Duhem inequality reads 

where p is mass density, and Q is the rate of entropy production per unit mass. Substi- 
tution of expression (22) into this equation implies that 



ToQ(t) = - [a{t) - piiA(t)] §(^) + ^ [^i(^) + ^(^) 



> 0. 



(26) 



Because Eq. (26) is to be fulfilled for an arbitrary program of straining, e = e(t), the 
expression in the first square brackets vanishes. This assertion together with Eq. (24) 
results in the stress-strain relation 



r r n 1 /"* on r n 1 

a{t) ^ E^[e{t) - e^{t)\ - - du | -^{t,T,u)[e{T) - e^{T)\dTj, 



(27) 



where E = ppX is an analog of the Young modulus. It follows from Eqs. (18), (23), (24) 
and (27) that 

= (28) 

According to Eqs. (25) and (28), the functions Ai(t) and A2{t) are non-negative for an 
arbitrary program of loading, which implies that the Clausius-Duhem inequality (26) is 
satisfied. 

Substitution of Eqs. (3), (14) and (15) into Eq. (27) results in the constitutive 
equation 



a{t) = e\ e{t) - eu(t) - «a(e(t)) / pico)du; 



X 



Jo 



exp 



r{e{s),uj)ds 



e(T) - eu(T) 



dry 



(29) 



Given functions p{w), ra(e) and n^i^), the time-dependent response of a semicrystalline 
polymer at isothermal uniaxial loading with small strains is determined by Eqs. (2), (18) 
and (29). For a standard relaxation test with the longitudinal strain e°, 



e{t) 



0, t<0, 
e°, t > 0, 
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these equations imply that 



{t, e°) = Ci(eO) - C2(e°) H p{u;) [l - exp(-r,(e°) exp{-u;)t)\dcu, (30) 



where e' 



is the strain induced by shding of junctions and 



C,{e') = E{e' - el) 



(31) 



To fit experimental data, we adopt the random energy model with 



p(uj) = po exp 



2S2 



uj>0 



p{u) = 0, a; < 0, 



(32) 



where and S are adjustable parameters, and the pre-factor po is determined by the 
condition _ 



Given a strain e°, Eqs. (30) and (32) are determined by 5 material constants: 

1. an analog of the average potential energy for rearrangement of strands Q, 

2. an analog of the standard deviation for distribution of potential energies S, 

3. the attempt rate for separation of strands from temporary junctions in active MRs 

4. the coefficients Ci and C2. 

Our aim is to determine these parameters by fitting experimental data depicted in Figures 
3 to 11. 

5 Fitting of observations 

We begin with matching relaxation curves for specimens not subjected to thermal tr- 
eratment. First, we approximate experimental data measured at the strain €2 = 0.04. 
This strain is chosen because it is located substantially below the yield point, on the one 
hand, and the testing machine ensures high accuracy of control of the strain level in the 
relaxation mode, on the other. 

Because the rate of rearrangement. Fa, and the average potential energy, Q, are mu- 
tually dependent [according to Eqs. (30) and (32), the growth of Q results in an increase 
in Fa], we set Fa = 1 s and approximate the relaxation curve by using 4 experimental 
constants: Q, S, Ci and C2. To find these quantities, we fix the intervals [0,f2max] and 
[0, Smax], where the "best-fit" parameters Q and S are assumed to be located, and divide 
these intervals into J subintervals by the points Qi = iA^ and Sj = j'As {i,j = l,..., J) 
with An = flraax/J, Aj] = Sjnax/-'^- For any pair, {Qi,T,j}, we evaluate the integral in 
Eq. (30) numerically (by Simpson's method with 200 points and the step = 0.1). The 




(33) 
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pre-factor po is determined by Eq. (33). The coefficients Ci = Ci{i,j) and C2 = €2(1, j) 
are found by the least-squares method from the condition of minimum of the function 



^-^(^) j) — ^ ^ '^exp(^m) '^num(^m) 



(34) 



'■m 



where the sum is calculated over all experimental points tm- The stress (Xcxp in Eq- (34) is 
measured in the relaxation test, whereas the stress (Xnum is given by Eq. (30). The "best- 
fit" parameters and S minimize the function J' on the set I Cli, (i, j = 1, . . . , J) >. 



After determining the "best-fit" values, fij and S^, we repeat this procedure for the 
new intervals and [Ej_i,Ej+i] to ensure good accuracy of fitting. Figure 4 

demonstrates fair agreement between the experimental data and the results of numerical 
simulation with Q = 4.29 and S = 4.34. 

To approximate relaxation curves at other strains, e^, we fix the constants Q and E 
found by matching observations at 62 and fit every relaxation curve by using 3 adjustable 
parameters: Fa, Ci and C2. These quantities are determined by using a procedure similar 
to that employed in the approximation of the relaxation curve at £2. We fix the interval 
[0,Fmax], where the "best-fit" attempt rate Fa is supposed to be located, and divide this 
interval into J subintcrvals by the points Fj = iA-p {i = 1, . . . ,J) with Ar = Fmax/^- For 
any Fj, we calculate the integral in Eq. (30) numerically and calculate the coefficients 
Ci = Ci{i) and C2 = C2{i) by the least-squares method from the condition of minimum for 
function (34). The "best-fit" attempt rate minimizes the function J on the set |Fj (i — 

1, . . . , J)|. When this "best- fit" value, Fj, is found, the procedure is repeated for the new 
interval [Fj_i, Fj+i] to ensure an acceptable accuracy of fitting. Figures 3 to 11 show good 

agreement between the observations and the results of numerical analysis. 

The above algorithm of fitting is repeated to approximate the relaxation curves for 
specimens annealed for 4 and 24 hours. The "best-fit" parameters fl and E read 5.70 
and 4.88 for samples annealed for 4 h and 5.19 and 3.80 for specimens annealed for 24 h, 
respectively. 

For a quasi-Gaussian distribution function (32), the parameters fl and E do not co- 
incide with the average potential energy for detachment of active strainds, Qq, and the 
standard deviation of potential energies for separation of strands from the network, Eq. 
The latter quantities read 



The dimensionless parameters Qq, Eq and ^ = Hq/Qq given by Eq. (35) are listed in 
Table 1 which shows that the width of the quasi-Gaussian distribution (characterized by 
the ratio ^) monotonically decreases with annealing time (however, changes in ^ are rather 
weak) . 

For any longitudinal strain e^, the attempt rate, Fa(efe), is determined by matching an 
appropriate relaxation curve. The fraction of active MRs, Ka(cfc)) is found from Eq. (31), 




(35) 



«a(efc) 



C2{ek) 



Ci(efe)- 
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These quantities are plotted versus strain e in Figures 12 and 13. The experimental data 
are approximated by the phenomenological equations 



logTa = 7o + 7ie, 



(36) 



where the coefficients 7^ and fcj are found by the least-squares method. Figures 12 and 13 
reveal that the viscoelastic behavior of specimens not subjected to thermal treatment is 
rheologically simple in the sense that the quantities Fa and Ka in Eq. (30) are independent 
of strain. On the contrary, annealed samples demonstrate the time-dependent behavior 
that is strongly affected by loading: with an increase in strain, the attempt rate. Fa, 
exponentially decreases and the fraction of active MRs, k^, linearly grows. 

According to Figure 12, the attempt rate. Fa, is strain-independent for non-annealed 



A detailed 
shows that the attempt rate 



specimens (however, the scatter of the experimental data is rather large), 
analysis of relaxation curves for non-annealed samples 
increases in the range of strains from e = 0.005 to e 
This implies that the free-volume concept 



m 



0.02 and remains constant at 
|, |U| is valid for isotactic 



e > 0.02. 

polypropylene, but the area of its applicability is confined to relatively small deformations 
far below the yield strain. 



6 Discussion 

Several approaches were recently proposed to the description of the viscoplastic behavior of 
isotactic polypropylene in the vicinity of the yield point. Aboulfaraj et al. |^ presumed 



that plastic slip mechanisms had noticeably different features in a and (3 structures. 
Karger-Kocsis and Varga [50| and Karger-Kocsis et al. |^| explained toughening of iPP 



by mechanically-induced /3 — > a transformation of crystallites. Raab et al. [q^ associated 
the difference in the response of a and (3 spherulites with different types of chain folding 
in lamellae. 

Three substantial shortcoming of these concepts should be mentioned: 

1. they are based on some hypotheses about the difference in the mechanical behavior 
of a and (3 crystallites which have not yet been confirmed experimentally, 

2. these models imply that changes in the stress-strain diagrams of iPP at annealing 
are associated with an increase in the content of /5-polymorph, which contradicts to 



WAXS (wide angle X-ray scattering) diagrams obtained by lijima and Strobl [|26 



3. these approaches do not expound a pronounsed decrease in the relaxation rate with 
strain for annealed samples exhibited in Figure 12. 

The results presented in Figures 1, 12 and 13 may be explained in terms of the growth of 
heterogeneity in the distribution of lamellar strength notwithstanding whether the fraction 
of (3 spherulites increases at annealing. 

An increase in the low-temperature shoulder of the DSC traces depicted in Figure 1 
(at a practically constant enthalpy of melting) means that the content of "weak" crys- 
talline lamellae (that melt at relatively low temperatures) noticeably grows. These "weak" 
lamellae may be associated with subsidiary lamellae in a spherulites. 
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Using a thermo-mechanical analogy, one can speculate that the growth of the fraction 
of thermally weak lamellae at annealing is tantamount to an increase in the concentration 
of lamellae which can be easily fragmented under stretching. The latter reflects an increase 
in structural disorder of a crystallites at annealing far below the melting point observed 
as substantial variations in the relative intensity of Bragg reflections . 



On the other hand, annealing of isotactic polypropylene induces thickening of radial 
lamellae, which results in the growth of elastic moduli of the polymer. This implies that at 
relatively small strains (e < 0.02 to 0.04, when fragmentation of lamellae does not occur), 
annealing of iPP leads to an increase in the longitudinal stress (which is demonstrated in 
Figure 2). At higher strains, the coarse slip starts in weak lamellae of specimens annealed 
for 24 h (which are less homogeneous than those annealed for 4 h), whereas no lamellar 
fragmentation takes place in samples annealed for 4 h. As a consequence, the stress-strain 
curve 2 is located higher that the curve 3 in Figure 2. With an increase in strain, lamellar 
fragmentation occurs in both specimens, which implies that the stress-strain curves for 
annealed specimens practically coincide at strains exceeding the yield point ey = 0.13. 

Further increase in strain in the post-yield region results in total fragmentation of 
lamellae that can be broken at a given stress intensity in specimens not subjected to 
thermal treatment. Because homogeneous (along a specimen) crystal slip becomes im- 
possible, spatial heterogeneity in the distribution of stresses arises at the macro-level, 
which leads to necking of the specimens at en = 0.18. For the samples annealed for 4 and 
24 hours, the concentration of "weak" lamellae is higher, which implies that the strains 
corresponding to the total fragmentation of weak lamellae and transition to the spa- 
tially heterogeneous deformation of specimens exceed that for the non-annealed material 
(en = 0.25 and en > 0.30, respectively). 

Noticeable fragmentation of lamellae in spherulites implies that some amorphous re- 
gions are released (whose deformation was previously screened by surrounding lamellae 
in non-broken crystallites), which results in an increase in the fraction of active MRs, Ka, 
with strain. This release of restricted amorphous phase is substantially less pronounced 
in the specimens not subjected to thermal treatment (curve 1 in Figure 13) compared 
with the annealed specimens (curves 2 and 3 in Figure 13). Three features of the curves 
depicted in Figure 13 are worth to be mentioned: 

1. At small strains (less than 0.04), the concentration of active MRs monotonically 
decreases with annealing time. This phenomenon may be attributed to an increase 
in the fraction of amorphous regions whose mobility is restricted by surrounding 
lamellae at annealing (driven by development of subsidiary lamellae). 

2. The content of active meso-domains linearly grows with strain in annealed specimens 
with the rate that is practically independent of the annealing time. This may be 
explained by changes in the micro-structure of spherulites at annealing: although 
the rate of lamellar fragmentation is assumed to be higher in iPP annealed for 24 h, 
the amount of amorphous phase released at any fragmentation act is smaller than 
in iPP annealed for 4 h. 



3. Curves 2 and 3 intersect curve 1 in the region between e = 0.08 and e = 0.12, i.e. 
in the close vicinity of the yield point for non-annealed specimens (see Figure 2). 
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With reference to Coulon et al. PS| and Raab et al. p2[ , we suppose that under stretching 



lamellae are fragmented into small aligned blocks that serve as extra physical cross-links 
in amorphous meso-domains. According to the concept of transient networks |^^, an 
increase in the concentration of permanent cross-links results in a decrease in the net 
rate of rearrangement. This rate is characterized in the model by the attempt rate, Fa, 
which is considered as an average (over MRs) rate of detachment of active strands from 
their junctions. This conclusion is fairly well confirmed by the results depicted in Figure 
12: the attempt rate is practically independent of strain for the specimens not subjected 
to thermal pre-treatment, and Fa exponentially decreases with strain for the annealed 
samples. Dispite apparent similarity in the slopes of curves 2 and 3 in Figure 12, it is 
rather difficult to assert that the kinetics of the strain-induced decrease in the attempt 
rate is independent of annealing time because of the large scatter of data for the specimens 
annealed for 24 h. It is worth noting that the attempt rates were determined in Section 
5 from the condition Fa = 1 s at e = 0.04, which implies that their values cannot be 
directly compared for specimens annealed for different amounts of time (because these 
samples have different distributions of potential energies for separation of active strands 
from temporary nodes, see Table 1). 



7 Concluding remarks 

Constitutive equations have been derived for the time- dependent behavior of semicrys- 
talline polymers at isothermal loading with small strains. A mean-field approach is 
employed to develop stress-strain relations: a complicated micro-structure of isotactic 
polypropylene is replaced by an equivalent transient network of macromolecules bridged 
by junctions (physical cross-links, entanglements and crystalline lamellae). The network 
is assumed to be strongly inhomogeneous, and it is thought of as an ensemble of meso- 
regions with various potential energies for separation of strands from temporary nodes. 

The viscoelastic response of a semicrystalline polymer is ascribed to separation of 
active strands from temporary junctions and merging of dangling strands to the network 
in active meso-domains. Rearrangement of strands is modelled as a thermo-mechanically 
activated process whose rate is given by the Eyring equation with a strain-dependent 
attempt rate. 

The viscoplastic response is described by slippage of junctions with respect to their 
positions in the bulk material. The rate of sliding is assumed to be proportional to the 
macro-stress in a specimen. 

Three series of tensile relaxation tests have been performed on isotactic polypropylene 
at ambient temperature. In the first series, injection-molded samples are used without 
thermal pre-treatment. In the second series, the samples are annealed at 130 °C for 4 
h, and in the last series, the speciments are annealed for 24 h at the same temperature. 
Adjustable parameters in the stress-strain relations are found by fitting observations in 
the range of strains from 0.02 to 0.18. The following conclusions are drawn from the 
analysis of experimental data: 

1. The relaxation spectrum of iPP (which is determined by the distribution function, 
p{uj), for potential energies of detachment of active strands from their junctions) is 
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not affected by mechanical factors, but is altered at annealing. 

2. The attempt rate, Fa, for separation of active strands from temporary nodes is 
practically independent of strain for specimens not subjected to thermal treatment, 
and it exponentially decreases with strain for annealed samples. 

3. The relaxation strength (which is characterized by the content of active mcso-regions 
Ka) is independent of strain for non-annealed specimens, and it linearly increases 
with strain for annealed samples. 

These findings are qualitatively explained based on the hypothesis that the distribution 
of strengths of crystalline lamellae in iPP is noticeably broadened at annealing, which 
results not only in thickening of lamellae, but also in the growth of "weak" (subsidiary) 
lamellae that are easily fragmented either by heating or by mechanical loading. DSC 
measurements provide some confirmation for this assumption. 
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List of figures 



Figure 1: DSC melting curves for a non-annealed specimen (unfilled circles), a specimen 
annealed for 4 h (filled circles) and a specimen annealed for 24 h (triangles). Symbols: 
experimental data 

Figure 2: The longitudinal stress a MPa versus strain e in tensile tests with the cross- 
head speed 5 mm/min. Symbols: experimental data. Unfilled circles: a virgin specimen; 
filled circles: a specimen annealed for 4 h; triangles: a specimen annealed for 24 h; 
asterisks: stresses at the beginning of relaxation tests at various strains 

Figure 3: The longitudinal stress a MPa versus time i s in a tensile relaxation test at 
e = 0.02. Symbols: experimental data. Solid lines: results of numerical simulation. Curve 
1: a virgin specimen; curve 2: a specimen annealed for 4 h; curve 3: a specimen annealed 
for 24 h 

Figure 4: The longitudinal stress a MPa versus time t s in a tensile relaxation test at 
e = 0.04. Symbols: experimental data. Solid lines: results of numerical simulation. Curve 
1: a virgin specimen; curve 2: a specimen annealed for 4 h; curve 3: a specimen annealed 
for 24 h 

Figure 5: The longitudinal stress a MPa versus time i s in a tensile relaxation test at 
e = 0.06. Symbols: experimental data. Solid lines: results of numerical simulation. Curve 
1: a virgin specimen; curve 2: a specimen annealed for 4 h; curve 3: a specimen annealed 
for 24 h 

Figure 6: The longitudinal stress a MPa versus time t s in a tensile relaxation test at 
e = 0.08. Symbols: experimental data. Solid lines: results of numerical simulation. Curve 
1: a virgin specimen; curve 2: a specimen annealed for 4 h; curve 3: a specimen annealed 
for 24 h 

Figure 7: The longitudinal stress a MPa versus time t s in a tensile relaxation test at 
e = 0.10. Symbols: experimental data. Solid lines: results of numerical simulation. Curve 
1: a virgin specimen; curve 2: a specimen annealed for 4 h; curve 3: a specimen annealed 
for 24 

Figure 8: The longitudinal stress a MPa versus time t s in a tensile relaxation test at 
e = 0.12. Symbols: experimental data. Solid lines: results of numerical simulation. Curve 
1: a virgin specimen; curve 2: a specimen annealed for 4 h; curve 3: a specimen annealed 
for 24 h 

Figure 9: The longitudinal stress a MPa versus time i s in a tensile relaxation test at 
e = 0.14. Symbols: experimental data. Solid lines: results of numerical simulation. Curve 
1: a virgin specimen; curve 2: a specimen annealed for 4 h; curve 3: a specimen annealed 
for 24 h 

Figure 10: The longitudinal stress a MPa versus time t s in a tensile relaxation test 
at e = 0.16. Symbols: experimental data. Solid lines: results of numerical simulation. 
Curve 1: a virgin specimen; curve 2: a specimen annealed for 4 h; curve 3: a specimen 
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annealed for 24 h 

Figure 11: The longitudinal stress o MPa versus time t s in a tensile relaxation test 
at e = 0.18. Symbols: experimental data. Solid lines: results of numerical simulation. 
Curve 2: a specimen annealed for 4 h; curve 3: a specimen annealed for 24 h 

Figure 12: The attempt rate Fa versus strain e in tensile relaxation tests. Symbols: 
treatment of observations. Unfilled circles: virgin specimens; filled circles: specimens 
annealed for 4 h; triangles: specimens annealed for 24 h. Solid lines: approximation of 
the experimental data by Eq. (36). Curve 1: 70 = —0.15, 71 = —0.08; curve 2: 70 = 0.12, 
71 = -2.91; curve 3: 70 = 0.05, 71 = -2.99 

Figure 13: The concentration of active MRs versus strain e in tensile relaxation tests. 
Symbols: treatment of observations. Unfilled circles: virgin specimens; filled circles: spec- 
imens annealed for 4 h; triangles: specimens annealed for 24 h. Solid lines: approximation 
of the experimental data by Eq. (36). Curve \: k^ — 0.57, kx — 0.01; curve 2: /cq = 0.47, 
kx = 1.55; curve 3: k^ = 0.39, kx = 1.49 
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Table 1: Adjustable parameters Qq, and ^ at various annealing times h 
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